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Abstract 

This paper is intended to investigate the relation between electrodynamics 
in anisotropic material media and its analogous formulation in an spacetime, with 
non-null Riemann curvature tensor. After discussing the electromagnetism via chiral 
differential forms, we point out the optical activity of a given material medium, closely 
related to topological spin, and the Faraday rotation, associated to topological tor- 
sion. Both quantities are defined in terms of the magnetic potential and the electric 
and magnetic fields and excitations. We revisit some properties of material media and 
the associated Green dyadics. Some related features of ferrite are also investigated. 
It is well-known that the constitutive tensor is essentially equivalent to the Riemann 
curvature tensor. In order to investigate the propagation of electromagnetic waves in 
material media, we prove that it is analogous to consider the electromagnetic wave 
propagation in the vacuum, but this time in a curved spacetime, which is obtained 
by a deformation of the Lorenztian metric of Minkowski spacetime. Spacetime de- 
formations leave invariant the form of Maxwell equations. Also, there exists a close 
relation between Maxwell equations in curved spacetime and in an anisotropic mate- 
rial medium, indicating that electromagnetism and spacetime properties are deeply 
related. For instance, the equations of holomorphy in Minkowski spacetime are es- 
sentially Maxwell equations in vacuum. Besides, the geometrical aspects of wave 
propagation can be described by an effective geometry which represents a modifi- 
cation of the Lorentzian metric of Minkowski spacetime, i.e., a kind of spacetime 
deformation. 
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Introduction 



The metric- free formulation of electromagnetism is an old concept E3 1281 
1311 1331 1551 110| . and it is more natural, correct, precise and geometrically sensible if 
differential forms, intrinsically endowed with chirality [161 I28| . are used |KSI I33| . The 
excitations 1 T)(x), H(x) (respectively the electric displacement and the magnetic field) 
and the fields ~E(x), B(x) (respectively the electric field and the magnetic induction) 
are naturally described in a spacetime destituted of a metric, which is introduced 
only when the constitutive relations are to be defined. The metric-free (and affine- 
free) formulation of electrodynamics brings a geometric character and a clear physical 
interpretation. But if we want to relate the Faraday excitation 2-form field G{x) and 
the electromagnetic field strenght 2-form field F(x), we need to consider a constitutive 
tensor \i that gives the relation G{x) = \F(x). Equivalently, spacetime must be 
endowed with a metric, since the constitutive law depends on the spacetime metric [5] 

The constitutive tensor (CT) is more than a relation between F(x) and G(x), 
describing physical intrinsic properties of matter or spacetime. It is essentially the 
Hodge star operator that, besides the well-known duality between fc-forms and (n— 
fc)-forms in a rt-dimensional vector space (endowed with a metric) , changes the parity 
of differential forms, but with an additional information about the medium structure. 
In this sense, CT describes the properties of spacetime (magnetic) permeability and 
(electric) permittivity, since a general medium can be arbitrarily anisotropic. Under 
this viewpoint, a CT immediately brings a light-cone conformal structure to spacetime 
|15| . Formally, the link between CT and spacetime structure is expected, since CT 
and the Rlemann curvature tensor have the same mathematical properties. A CT can 
reveal precious informations about spacetime, for example, the CT scalar curvature is 
identically null in any medium possessing central symmetry |33| . 

In this paper the CT \ that describes any linear (in particular, crystalline) media, 
in the general case presenting optical activity, is expressed by a conformal transfor- 
mation of the vacuum CT xo- The metric associated to the medium is derived from 
the CT (up to a conformal factor) and vice- versa |2(J1 1221 |2"31 125) . In this sense, the 
Lorentzian metric of Minkowski spacetime, associated with xo: is deformed into a 
general metric of a riemannian spacetime, related to x- 

We prove that in order to completely describe the CT of any linear medium pre- 
senting natural optical activity, we only need the matrix 7, that describes the optical 
activity of such a medium, and xo- ln the particular investigation of crystalline media, 
we describe the constitutive tensor associated to the 32 crystal classes presenting nat- 
ural optical activity uniquely from xoj i-e., from the spacetime metric, since xo can be 
written as second order metric tensor combination |22l I2(J| . This paper is organized 
as follows: in Sec. 1 we review the geometric description of fields and excitations 
J We follow the nomenclature given in \22\. 



in electromagnetic theory, well-established originally in the papers by Maxwell \21\ . 
Schonberg QU], Hehl [22], Kiehn gU, Post gSj and Jancewicz J2HI2B], among others. 
The Maxwell equations carry information about the nature of the different fields and 
excitations, and the theory arises with an essential geometric character, if the fields 
and excitations are correctly interpreted, as in |27l 1281 \T2\ as even and odd dif- 
ferential form fields. In Sec. 2 we revisit the main features of electromagnetism in 
anisotropic material media. The Green dyadic is obtained and we treat the particular 
case of a electrically anisotropic material media. Plane waves are investigated, illus- 
trating the present approach. Also, optically active media are treated. For instance, 
ferrite is explicitly studied. In Sec. 3 Maxwell equations in curved spacetime and in 
an anisotropic material medium are proved to be equivalent. It sheds some new light 
on the differential geometric aspects of electromagnetic wave propagation, that can be 
described by a geometry which represents a modification of the metric in Minkowski 
spacetime. We prove that all the information contained in the constitutive tensor as- 
sociated to any linear media (arbitrarily describing optical activity), is precisely given 
by the vacuum constitutive tensor and by the matrix that describes the optical activity 
associated with the medium. Such matrix is essentially viewed as a perturbation of 
the original constitutive tensor that does not describe optically active media. Finally 
in the Appendix, the main results concerning differential forms are reviewed. 

1 Electromagnetism in the exterior algebra 

Heretofore Vt k {M) denotes the space of fc-form fields defined on a manifold M. 
Given the map E : M — > R 3 , the electric field E(x), x € M, is an even 1-form 
field (E(x) S i7 1 (]R 3 )), since E(x) is a linear map from the infinitesimal vector dr(x) 
to the infinitesimal scalar potential dV(x), given by dV(x) = — E(x) • dr(x). The 
physical dimension of E(x) in the SI, [E(x)] = Vm _1 , agrees with this interpretation. 
Analogously the magnetic induction B(x) is an even 2-form field (B(x) € f2 2 (R 3 )), 
since B(x) is a linear map from the infinitesimal bivector dS(x) to the infinitesimal 
scalar d<p(x). Explicitly we have d<f>(x) = — B(x) • dS(x), where <j)(x) is the magnetic 
flux. The physical dimension B(x) in the SI, [B(x)] = Wbm -2 = T ( = Tesla), again 
agrees with such an intepretation |21l |2Hj • 

From now on we call an even (odd) differential form field the one that doesn't (does) 
change sign under parity transformations 2 |16l I28j . Even form fields are elements of 
£2+ (M) , hereon simply denoted by fi(M) , while odd form fields are elements of f2_ (M) . 
Such forms are called chiral differential forms. 
2 A parity transformation is defined in R n as the inversion of an odd number of basis vectors of K n . 



1.1 Homogeneous Maxwell equations and potentials 

The differential operator d : Q±(M) — > 05_ +1 (M) does not change the differential 
forms parity. From now on we denote d a = d/da. We also adopt natural units, in 
particular the speed of light c = f , in what follows. 

The first homogeneous Maxwell equation is an expression relating even 2-form 
fields: 

dE(x) + d t B{x) = 0. (1) 
The absence of magnetic monopoles can be described by the equation 

<m(x) = o. (2) 

Eqs.JU and J2J are the homogeneous Maxwell equations. From eq.J3J), using the 
Poincare lemma 3 , there exists an even 1-form field A(x) (the magnetic potential) 
satisfying the relation 

B(x) = dA(x). (3) 

Substituting in eq.Q, one obtains the expression dE(x) + dtdA(x) — 0, or d (E(.x) + 
dtA(x)) — 0. Using the Poincare lemma, there exists a scalar field potential <&(:r) G 
n°(R 3 ) such that E(x) + d t A(x) = -d$(x), implying that 

E(x) = -d t A(x) -d$(x). (4) 



1.2 Non-homogeneous Maxwell equations 

The electric current density j(x) is an element of ^(M 3 ), an odd 2-form field, 
which changes sign under parity transformations. It is clear that dj e f7 3 _(IR 3 ), and 
then the continuity equation (local form of electric charge conservation) can be written 
as 

dj(x) + d t p(x) - 0, (5) 

where p(x), the electric charge density, is an odd 3-form field. Obviously dp(x) = 0, 
and the Poincare lemma again asserts that there exists D(x) S ^^(R 3 ) such that 

p{x) = d~D{x). (6) 

The 1-form field T>(x) is called electric displacement. The unity of D(x) in the SI is 
\D(x)] = Cm -2 (C = Coulomb). Substituting in eq.©, we have dj(x)+d t (dT>(x)) = 0, 
and d(j(x) + d t T)(x)) = 0. There exists an odd 1-form field U(x) E ^^(K 3 ) such that 



j(x)+d t T>{x) = dH(x). (7) 

3 The Poincare lemma asserts that if an open set 13 C R n is star-shaped, every closed form is exact in 
15. (uj £ Q, k (M) is closed if div — 0, and exact if there exists a form tj £ such that u) = drj. An open 

set 15 C R n is star-shaped with respect to the origin if, for all x € 13, the line from the origin to x is in 13.) 



It describes the Ampere- Oersted law. The SI unit of H(x) is [H(x)] = Am -1 (A = 
Ampere). The odd form fields D(x) and H(x) are potentials with sources p{x) and 
j(x), respectively. 

The Poyinting vector S(x) describes the electromagnetic strenght energy flux den- 
sity. It is possible to write 

S(x) =E(x) AH*(x). (8) 

From the algebraic viewpoint, the product above is the unique possibility, since 
quantities representing flux densities are described by odd differential 3-form fields 
1221 ESI HID EDO Ei2|, elements of ft 3 (R 3 ). The electric (w e (x)) and magnetic (w d (x)) 
energy densities are elements of fi^R 3 ), expressed by 

w e (x) = ^E(x)AD(x) = ~D(x)AE(x), w m = ~B(x)AH(x) = iH(x)AB(x). (9) 

The electromagnetic field energy density is written as [2S| 

w = w e + w m = -(E(x) A D(x) + B(x) A H(x)). (10) 

1.3 The Hodge star operator 

It is well-known j2Hl that the vector spaces fi fe (R 3 ) and f2 3 ~ fc (R 3 ) have the same 
dimension, since dim f2 fc (R 3 ) = ( 3 ) = ( 3 _ fc ) = dim £l 3_fc (R 3 ). The same result is 
valid to any n-dimensional space (see Appendix). Meanwhile, it docs not exist any 
canonical isomorphism between f2 fe (R 3 ) and S! 3 ~ fc (R 3 ). The isomorphism given by the 
Hodge star operator* : K k (V) — > A n ~ k (V) always satisfies ** = ±id. For more details, 
see the Appendix. The contraction is a generalization of the interior product, and it 
can be written in terms of the Hodge star operator and the exterior product, as 

i>j<P = *((*</>) A i>), (11) 

where ip,cj) € il(M 3 ). For more details see, e.g., [361 I37| . 

1.4 The Poynting theorem 

If we take the C-conjugation of eqs.JQ) and Q and respectively multiply by E*(x) 
and H*(x) we obtain 

E(x) A dH*(x) - H*(x) A dE(x) = - |e(x) A (d t *eE*(x)) + *u _1 B(a;) A (d t B*(x)) 

c L 

= *[e(x)j (9 t eE*(x)) + (/i" 1 B)(a;)j (d t B*(x))](.12) 
Ea. (|12|l can be written in a coordinate system as 

i<9 t ^E(x) j(eE(x)) + B(x) j(/i~ x B(x))) dx A dy A dz. 
and from the expression 



-dS(x) = H*(x) A dE(x) - E(x) A dH*(x), 



(13) 



it follows that 

dS(x) — dt'w(x), 
the so-called Poynting theorem |2l ll3llTl] . 



(14) 



1.5 Electromagnetic Intensity and Excitation 

The electromagnetic field strenght F(x) € f2 2 (R 1,3 ) is an even 2-form in R 1 ' 3 , also 
called the Faraday 2-form field 19 . If an arbitrary, but fixed, time vector is chosen 
in R 1,3 , we can split spacetime in space plus time. Then it is possible to use E(x) and 
B(x) to describe F(x) as 

F(x) = B(x) + E(x) Adt. (15) 

The electromagnetic excitation G{x) € fi?_(R 1,3 ) can also be considered as an odd 
2-form field given bv 1251133] 

G(x) = D(x) - H(x) A dt. (16) 
Eas. H1I2|I can be summarized as 

dF{x) = 0, (17) 

and eas. 1)617(1 are synthetically written as 

dG{x) = J(x), (18) 

when the odd 3-form current density field J{x) = p(x) — j(x) A dt is defined [2E]- If 
we admit primarily ea. l(17|l . the electric and magnetic fields are only defined after a 
spacetime splitting. 

1.6 Vacuum constitutive relations 

Hereon it is assumed implicitly that the Hodge star operator changes the parity 
of the differential forms 4 . 

Constitutive relations are written as 

G(x)=*F(x). (19) 

This relation can be expressed in the vacuum, after a spacetime splitting, as: 

D(x) =e Q *E(x), B(x)= /j,o*U(x), (20) 

where £o denotes the vacuum electric permittivity and /io denotes the vacuum mag- 
netic permeability. From ea. (|17fl it is possible to find A(x) G 57 1 (R 1,3 ) such that 

F(x) = dA{x). (21) 
4 This Hodge star operator is, de facto, the composition of the Hodge star operator with a pseudoscalar 
|42l 1431 This new operator is then able to lead odd (even) form fields to even (odd) ones. (By abuse 
of notation we also denote this new operator by *. 



The even 1-form field A(x) denotes the well-known electromagnetic potential. In 
components, ea. (|2~T]> is written as 

F^^dnAv-dvAp. (22) 

Ea. (|21|l is invariant under the maps A(x) i— > A(x) +u>(x), u>(x) G 51 1 (M 1 ' 3 ) such that 
du(x) — 0. In particular ea. (j21j) is invariant when A(x) i— ► A(x) + d\(x), A(x) G 
1 (M 1 ' 3 ). The existence of form fields that are closed, but not exact, gives rise to the 
physical monopole and solitons in fluids, concerning paramount and striking appli- 
cations such as superconductivity, topological defects and turbulent non-equilibrium 
thermodynamics of fluids, exhaustively investigated by Kiehn |31|. 

One [3T] defines the odd 3-form field topological spin S(x) = A(x) A G(x) G 
fi^QR 1 ' 3 ) and the even 3-form field topological torsion T(x) = A(x)AF(x) € ^ 3 (R^ 3 ). 
It can be shown that optical activity is closely related to topological spin, while Fara- 
day rotation is associated to topological torsion |31| . 

Under a spacetime splitting it can be seen that 

T(x) = A(x) A F(x) = (A(x) - 4>(x)dt) A (B(i) + E(x) A dt) 

= A(x) A B(x) + (A(x) A E(x) - (f>(x)B(x))dt (23) 

and 

S(x) = A{x) A G(x) = (A(x) - 4>{x)dt) A (D(a:) - H(x) A dt) 

= A(x) A D(sb) + (A(x) A H(x) - <f>(x)T>(x))dt (24) 

Kiehn [3T] shows that T(x) is related to the helicity, while S(x) is associated to 
chirality of the electromagnetic fields. The 3-form field energy-momentum, is defined 
if an arbitrary direction ej is chosen: 

U t (x) = i[F(x) A (e.^G(x)) - G{x) A (e^Ffr))]. (25) 

The 3-form field energy- momentum is invariant under pseudodual maps F(x) i— > 
(p(x)G(x) and G(x) i— » —F(x)/tp(x), where y(a;) is an arbitrary scalar field non-null 
in all points of M 1 ' 3 . 



2 Revisiting electromagnetism in non-homogeneous 
media 

In the last decade, a lot of manuscripts have been concerning electrodynamics in 
material media via differential forms. For instance, see E] 03 [7\ |37| ■ 

For permeability and permittivity tensors such that the product e fj, is diagonal- 
izable, the expression for the Green diadic EJ E| is given by 

/ exp(imif) \ 



dct ft 
Airr 



\ 



exp(im2r) 

expiimsf) j 



(26) 



( 2 2 2 \ 1 

+ + j^P, and mi, m 2 , TO3 denote eigenvalues associated 

with the matrix £*/^ _t , such that Re mi > 0, i=l,2,3. /ii, ^2, A*3 denote eigenvalues 

associated with /x. 

The diadic given by ea.(|26|l can be immediately written as the tensor product @JE]: 
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2.1 Material media 



^ f e lmif dx 1 dx 2 + e m ^dy x ® dy 2 + e^dz 1 $ dz 2 ") . 
r \ / 



(27) 



A medium is completely isotropic if the electric permittivity e and the mag- 
netic permeability [i can be written, respectively, as fi = diag(//i, /xi, fix) and e = 
diag(ei, £1, £1). If the elements of fi and £ in an anisotropic medium are hermitian, 
such medium is called eletrically or magnetically girotropic. For instance, a plasma 
with static magnetic field in the z axis 



(28) 
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is eletrically girotropic. A magnetic girotropic medium is given by 
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(29) 



2.2 Electrically anisotropic media 

Suppose that, in a given material medium, e — diag(£i, £2, £3) and fj, = I, where / 
denotes the identity matrix. It is immediate that f = far and det fi = /if, from where 
it can be seen that £/i~* = fa 1 diag(£i, £2, £3). Since the Green form g is given by @] 



g = [(det - w 2 £At _t ] 



then the diagonal components of g are given by 

9jj(k) - 



A solution of eq. ijSQt is given by 
/ 

Mi 

9 
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(32) 



exp(iw v /£i7iTr) J 



and then the Green diadic is given by 

g(ri, r 2 ) = -— ( exp(iuj^/einir) dx\ ® cfe 2 + exp(iw-v/£ 2 ^ir) dyi ® <i?/ 2 
47rr V 

+ exp(ia7^/£3/iir) dzi ® dz 2 ^j (33) 

where r = |jr*i — r 2 ||. This equation is originally obtained by Warnick [3]. 

In the particular case of an uniaxial medium, with e\ = £ 2 = e and £3 = e z , ca. (1331 
is led to 

ffffi, r 2 ) = - — ( exp(iw-^/£/iir) ctei ® <ix 2 + exp(ia;-i/£/iir) <iyi ® rfy 2 
47rr V 

+ exp(iw-y/£ z /iir) <fei ® dz 2 ^ (34) 

In any material isotropic medium we have e z = e, and the vacuum is obtained when 
£ 1— > £0 and /x 1— * /.to- In this case, 

< 35) 

2.3 Plane waves 

Heretofore we denote E = E(x), B = B(x), . . . , in order to simplify the notation to 
be used. Suppose that a plane wave propagates in the s direction and let the electric 
field be expressed by 

E = E exp(i(fcs • r — tot)). (36) 
Denoting n = k/u> = 1/v, eqs.Q and J7J gives 

D = — nH A s, B = nE A s. (37) 

Using the constitutive relation B = *v,H, eas.l|37[l give 

o 

71 

D = — *(sAE)As 
A 4 

2 r 

E-s(s-E) . (38) 



— * 



By abuse of notation, here /it denotes the magnetic permeability, a real number, instead 
of the tensor /i. Define the component Ej_ of E such that Ej_ • s = 0. Then, since Ej^ 
is in the plaquette defined by E A s, ea. (|38[) is written as 

D = — Ej_. (39) 
A* 

2.4 Presnel equations 

Now let £1, £ 2 and £3 be the eigenvalues of the matrix representation of e. From 
ea.i|38fl we have 

Ei = ^^E • a, (40) 



from where we obtain 



-E • s ■ 



-E • s 



n z - fie 3 



-E • s = — E • s 



3 2 
87 



E 



Since s ■ a = 1, then X^=i s 2 = 1, and from ea. i|42|) it follows that 



3 2 



i=l n s fid 



(41) 



(42) 



(43) 



If we define the so-called principal propagation velocity := (/X£j 1 ^ 2 , 0, 0), ea. (|43fl 
is lead to 

D 2 2 2 

(44) 



+ 



9 9 9 9 9 9 
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= 0. 



Eqs.|@3), and 1^3} are called Fresnel wave equations [T4] , 



2.5 Ferrite 

Ferrite is a material medium defined by 

/ a -i/3 (A 



/i = Mo 



— : /Lto/ir 



(45) 



i/3 a 
v 7/ 

From now on we consider j3 < a, and it is easy to see that det /i = ^(a 2 — /3 2 )7 and 
that the eigenvalues of [i r are a + (3, a — (3 and 7. From the expression 



(a 2 -/? 2 )2 7 2 Mo ( 



j/ f_\ 



\/a + (3 y/a- f3 yfi) 



it follows that 



/? 2 )M(- 
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+ /3 a- 13 7 



(46) 



The square roots of the eigenvalues of ui Sfi are given by 



mi = lu 



7712 = LO 



7 m — , m 3 =Lj 



— . ( 47 ) 



(a + (3) no ' 

where e is the diagonal element of the ferrite permittivity tensor, g is given by 
u 2 (a 2 -/? 2 )V 



Aur' 



where f = ^/ (a 2 — /3 2 )~/fiQr' and r ' = 
the expression 



75T?> 7^?' 7?)- Eq.63 is equivalent to 



(49) 



where go = ^{a 2 — f3 2 )^j^ /Airr' , r' = \\fl ' — ri '|| and r( are analogously defined as 
r'. When (3 = and a — 7 in ea. 148(1 it follows that 



2 2 

.9 = -: exp(iw^e/i ar)7, 

47rr 



where = ^/^i + + z ^ ■ Denoting /i = a/^o ea. (|50(l can be written as 

g = - — exp(iuiy/ejlr)l, 
which is the well-known expression for an isotropic medium. 



(50) 



(51) 



2.6 Faraday rotations 

From eq.J7|) and eci. (|36[l . it follows that 

-kdH = —iujeE. 
Taking the differential of the last equation we obtain 

d * dH = uj 2 e *h H 



Now, if we solve eas. <|37|) for B, we obtain 



B — */, , 



(52) 



(53) 



(54) 



where = H — (s • H)s. In components, the field given by eq.JMJ) is written as 



Bi = 



n 2 — c 2 £/x. 



■(s-H)si. 



(55) 



From ea. l|55(l the component of H in the z-direction is zero, and if we make the 
assumption that H = (Hidx + H 2 dy)e t( * kz ~ ult \ ea. f53"| gives 



a —i(3 



I ! 1 ' \ 9 / H\ 

1 1 _ y. 1 



\i(3 a J \H 2 J 

which solution is given by 

k 2 ^ = w 2 e/io(a + /3), k 2 _ = u 2 efio(a — /?), 



(56) 



(57) 



describing two (left- and right-handed) circularly polarized plane waves. Now, substi- 
tuting eqs.JSZI) in eq.J23) it follows that Hi = ±iA, if H 2 — A, A e C. The general 
solution of the system is 



H 



i Cl Ae ik + z + ic 2 Ae ik ~ 



dAe ik - z + c 2 Ae 



ik—Z 



Choosing c\ — —c 2 = |, we obtain 
w 2 2 



H 2 (z) = de lk + Z -i^e lk ~ z , 



.A 

< 

2 



l dy. (58) 



(59) 



which can be written as 

'k, - k 



H x (z) = A cos ^ — — — zj exp(-i(fc + + fc_)z/2), 
ff 2 (z) = A sin( fc+ ~ exp(-t(fc + + fc_)z/2). (60) 



Let SeR such that 

tan0=^=tan(^— (61) 

It is immediate that 

/ fc, — fc_ \ 

0fe = ( z J + 2fc-7r, fc is an integer. (62) 

Restricting 8 £ [0, 27r) it is clear that the phase difference between the left- and right- 
handed components is 29q, where 



r/„ = ^-zuj^/e^a + — — \ . (63) 



Consider k± in ea. l|57|l given by a second-order approximation, i.e., 

a la 2 ( a 
2/3 + 8/32 I/? 



k ± =uje^\l± — + -— + 0(-) ). (64) 



Substituting in ea. i|63|) we have 

6> = -^zuj^/a^e. (65) 
It shows the well-known result asserting that ferrite is indeed a non-reciprocal medium. 

3 Material media viewed as spacetime deformations 
in vacuum 

In the formalism in, e.g. |2S1 EH that describes the electromagnetism in linear 
media, the dual Hodge operator action is equivalent to the constitutive \ tensor action 
on 2-form fields: 

*A = xA, A £ O^l 1 ' 3 ). (66) 
If cartesian coordinates are introduced, ea. H66|) is equivalent to 



*A = \^ vaKX aKTp A^dz" A dx v . (67) 



3.1 The constitutive tensor 

In this subsection we present and discuss the main results in, e.g., HHHES], con- 
cerning the relation between the constitutive and the Ricmann curvature tensors. In 



linear media, the 2-form electromagnetic intensity F(x) € r2 2 (K 1,3 ) is related to the 
electromagnetic excitation G S fi^R 1 ' 3 ) by the equation 

F(x) = X G{x). (68) 

Using cartesian coordinates, F{x) and G(x) are expressed as 



G{x) = ^G^^dx** A dx v , F{x) = ^F,, u (x)dx^ A dx v (69) 



and F llv (x) and G )JiV (x) are related by 

G^(x) = \e^ aP X apaX F aX (x), (70) 
where e is the Levi-Civita tensor. The symmetry 

y\v<JK y\VK(J j^XvtTK yl>\<7K (7l) 

arises, since F^^x) and G^ v (x) are antisymmetric. Besides, the lagrangian density 
C{x) = G{x) A F(x) is written as 1531154] 

C(x) = \x X '" TR Fxu{x)F ITK {x). (72) 



From the relation 



there exists the relation 



9F Al ,(x)aF ffK (a;) flF M (s)fl^(x) 
= x aKXv (74) 

For more details, see |34| . The number of independent coordinates is 21, (using 
eas. H7H74f) '). which comes from the analogy to the Riemann curvature tensor. Only 
the antisymmetric combinations are non-trivial. An order two antisymmetric tensor 
has exactly six components and then there would exist 36 componentes. Expressing 
X as a 6 x 6 matrix (x € Hom(R 6 ,R 6 )) and using the bivectorial notation, (where the 
indices /, J, . . . = 01,02,03,23,31,12 are defined) we can see from ea. (l74H that the 
matrix x IJ is symmetric (x IJ — X JI ), an d there exists (6x7)/2 independent compo- 
nentes in x- Using physical arguments, one can show that in uniform media we have 
the relation [531134] 

X [XvaK] = 0. (75) 

In vacuum x can be written as 

X XvaK = Yoyfg{g x °g™ - g XK g v °), (76) 

where g is the determinant of g^ u and Y" is the vacuum admittance. 
From the constitutive tensor x Post defines two invariants: 

Xi = Xx Xv , (77) 



which is called scalar curvature of the medium described by x, and 



X2 = ewX^X^W*- (78) 

Post [Po72] proves that \i is non-zero for any medium and xi = x' AiATfi ' is identically 
null in any medium possessing central symmetry. The constitutive tensor is explicitly 
represented by: 
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The matrix \iik is the magnetic permeability matrix, sik is the electric permittivity 
matrix and 7^ is a matrix that describes the electric and magnetic polarization effects. 
One can prove that in media possesing central symmetry, the matrix 7^ is null 331 I34| . 
In isotropic media the relations 

7ifc = 0, eik = s Q 6 lk , Q k = 1 6ik, (79) 

are satisfied. In this case, \i = an d X2 = —^-^0/^0 [HH!< 

We shall study the light propagation in crystalline media presenting optical activity, 
which are characterized by 32 classes Each class is represented by a symmetry 

represented in the table: 
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2:2, 2/2,5 



C denotes central symmetry, S is the cyclic permutation of the indices, E x is a reflec- 
tion with respect to the yz plane (analogous definition for E y and E z ), S x is a rotation 
using the a; axis, followed by a reflection related to the yz plane (And analogous defi- 
nitions for S y and S z ). 

All crystal classes described in the above table present natural optical activity, 
and the corresponding respective matrices "/ki, composing the tensor x, are described 



below (the number before the matrices indicates the class number above described): 
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The matrices corresponding to the classes 29,32 describe the chiral vacuum |31j . 



4 Arbitrary constitutive tensors from the vacuum 
CT 

From now on we adopt the notation F = F(x), G = G(x), . . ., omitting the argu- 
ment x. 

► Spectral Theorem: There always exists a conformal transformation that di- 
agonalizes the constitutive tensor x ^ 

Considering the splitting US 1,3 ~ R 3 x R, we can write 

where e = So I and fi^ 1 = Hq 1 !- By the theorem above, there exists a matrix F 
composed by the eigenvectors of \ such that r _1 %r = A is a diagonal matrix. 

4.1 The chiral vacuum 

In order to illustrate the general approach, we firstly consider the chiral vacuum, 
described by the matrix 



-£ 70 



(81) 



6x6 

where 7^ = 711 /. The matrix l|81|) has eigenvalues 01, cr 2 and eigenvectors {(0, 0, <j\, 0, 0, 1), 
(0, o-i,0, 0, 1, 0), ((71, 0, 0, 1, 0, 0), (0, 0, o 2 , 0, 0, 1), (0, o 2 , 0, 0, 1, 0), (o 2 , 0, 0, 1, 0, 0)}, where 



^1,2= Mo 1 -e ±^ . (82) 

2711 



Then 

r- 1 Y r = A = I 

s 2 



r i Ar = A=(- Sl ] (83) 

6x6 

where Si = —<j\ I and £2 = <?2 /. Denoting F = and G = T~ 1 G, we obtain 

(84) 

Define the odd form fields D and H as 

D = E^D, H = E^H. (85) 



It follows that 



(86) 



and we prove that 



G = XoF, (87) 

where xo vacuum constitutive tensor 5 . We obtain for the chiral vacuum, after doing 
the inverse maps, the constitutive relation 

G = X °F (88) 

where 

X ° = TAxor- 1 (89) 

Then the constitutive tensor x'^ ; related to the chiral vacuum, is completely described 
by the matrix 7 and the vacuum constitutive tensor. We only used conformal maps in 
R 1 ' , which are elements of the group. This kind of structure in electromagnetism was 
discovered by Bateman [35], who was the first to observe that the Maxwell equations 
are invariant under the conformal group [311 132j . 

4.2 Arbitrary linear media: crystalline media, optical activity, 
magnetic and dieletric Faraday effects 

The method is fundamentally analogous to the chiral vacuum case. Consider an 
arbitrary linear media described by the matrix 

X=(~" l) , (90) 



3 modulo dilation of the axis e^, es and e6 by /xo and contraction of ei, e2 and e3 by eo- 



/- 



A 



where 

' 711 712 7i3 

7 = 721 722 723 
y731 732 733 J 

The matrix H9U|) has eigenvalues a a [A = 1, 2, . . . , 6). Then 



(91) 



r- 1 x r = A 

where Si = — diag(<7i, a%, a 3) and £2 = 
G = r _1 G, we obtain 



-Ei 
S 2 , 



(92) 



-diag((T4, 05, ctq). Denoting F = T 1 F and 



E 
B 



-E 
B 



= r-^r 



Defining the vectors D = T, 1 1 D, H = S 2 1 H, it follows that 



' 



-D 
H 



-1 0^ 
1/ 



E 
B 



-E 
B 



-Si 
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(93) 



(94) 



This implies that 

G = XoF, (95) 

where \a is the vacuum constitutive relations. Calculating the inverse maps, we obtain 
for any medium the constitutive relation 



where 



G = X F 



X = TAxoT- 1 



(96) 



(97) 



The constitutive tensor associated to the 32 crystal classes presenting natural optical 
activity is described uniquely from xoi i-e., from the spacetime metric, since 



\V(JK 

AO 



Using coordinates we write 



(98) 



Note that the expression above is the constitutive relation for any crystalline material, 
and it depends only of the matrix 7 (given at the end of Sec. 3, for all crystal classes), 
that describes optical natural activity. Then it can be seen as the deformation of the 
metric of Minkowski spacetime into a metric of curved riemannian spacetime, since 
in order to describe the constitutive relations of any crystalline medium we only need 
the metric of Minkowski spacetime. 

In particular, it is also possible to express, from the Lorentzian metric of Minkowski 
spacetime, the constitutive tensor associated to the dielectric and magnetic Faraday 



rotations, and the natural optical activity in arbitrary rotational symmetric media. It 
is respectively given by the following matrices: 
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(100) 



where £ij = (^ 1 )y. Post proves [Po97] that electromagnetic waves propagate with 
phase velocity u given by 

u = ±(( £ ±£ 23 vr l/2 , (ioi) 

(dielectric Faraday rotation), 

u = ±^^, (102) 

(magnetic Faraday rotation) and 

« = ±^±J^ + 4. (103) 

£ y £/X £^ 

(natural optical activity). 

In the whole process described in this subsection, we only have accomplished con- 
formal transformations in R 1 ' 3 . 

Concluding Remarks 

We investigated the relation between electrodynamics in anisotropic material me- 
dia and its analogous formulation in an spacetime, with non-null Ricmann curva- 
ture tensor. The propagation of electromagnetic waves in material media is proved 
to be analogous to consider the electromagnetic wave propagation in the vacuum, 
now in a curved spacetime, which is obtained by a deformation of the Lorcnztian 
metric of Minkowski spacetime. Such process of performing deformations of the 
metric of Minkowski spacetime can be rigorously described using extensors. Also, 
there exists a close relation between Maxwell equations in curved spacetime and in an 
anisotropic material medium, indicating that electromagnctism and spacetime prop- 
erties are deeply related. Besides, the geometrical aspects of wave propagation can be 



described by an effective geometry which represents a modification of the Lorentzian 
metric of Minkowski spacetime. 

We discussed the optical activity of a given material medium, closely related to 
topological spin, and the Faraday rotation, associated to topological torsion. Both 
quantities are defined in terms of the magnetic potential and the electric and mag- 
netic fields and excitations. The existence of form fields that are closed, but not exact, 
gives rise to the monopole and solitons in fluids, concerning topological defects and 
turbulent non-equilibrium thermodynamics, exhaustively investigated by Kiehn 
In a forthcoming paper, since the integral over K 3 of the topological torsion spatial 
component, introduced by eci. (|23[) . is the writhe of a framed oriented link, it is possible 
to investigate link invariants in gauge theory, from the knot theory viewpoint. 
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